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Abstract 

In the problem of motion of the Kowalevski top in a double force field the 4-dimen- 
sional invariant submanifold of the phase space was pointed out by M.P. Kharlamov 
(Mekh. Tverd. Tela, 32, 2002). We show that the equations of motion on this 
manifold can be separated by the appropriate change of variables, the new variables 
si,S2 being elliptic functions of time. The natural phase variables (components of 
the angular velocity and the direction vectors of the forces with respect to the 
movable basis) are expressed via si,S2 explicitly in elementary algebraic functions. 
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1 Introduction 



The famous solution of S. Kowalevski pQ for the motion of a heavy rigid body 
about a fixed point was generalized for the case of double constant force field 
in [2|3] . This Hamiltonian system essentially has three degrees of freedom, and 
hardly can receive a clear geometrical or mechanical description of all types 
of motions. Invariant subsystems, which can be interpreted as systems with 
two degrees of freedom, were found in [2~f4"] . Phase topology of the case [2] was 
studied in [5]. The present paper deals with the case [I]. We show that by 
proper choice of local coordinates it is possible to obtain separated differential 
equations of motion, and express all phase variables explicitly in terms of two 
new variables, the latter being elliptic functions of time. 
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2 Equations of motion and known first integrals 



Consider a heavy magnetized rigid body with a fixed point placed in gravita- 
tional and magnetic constant force fields. Let a, (3 be the direction vectors of 
the force fields and e*i,e*2 be the radius vector of mass center and the vector 
of magnetic moment of the body. The scalar characteristics (for example, the 
product of weight and distance from the mass center to the fixed point) may 
be included in the length of either vector of the associated pair. We prefer to 
consider e\,e 2 to be unit vectors according to the model accepted in [3]. The 
Euler-Poinsot equations of motion have the form 

r dul r ^ ^ ^ _ z da _ d(3 -> 

I — = Iuj x uj + ei x a + e 2 x (3, — = a x to, — = (3 x u, (1) 
dt dt dt 

where uj is the angular velocity, I is the inertia tensor. All vector or tensor 
objects are referred to the basis attached to the body. 

Suppose that the body is dynamically symmetric and denote by 7r e the equa- 
torial plane of inertia ellipsoid. Choosing measure units one can always make 
the moment of inertia with respect to symmetry axis equal 1. Let 



J = diag{2,2,l}, 

ei,e 2 6 7r e , e x ■ e 2 = 0. 



Then ei, e*2 may be taken as the first vectors of the movable basis. It is known 
that under conditions (2) the system (1) is completely integrable due to the 
existence of the first integrals [2] 

H = u 2 l +ul + \ul-{a 1 + h) ) 
K = (uf -ul + ax- (3 2 ) 2 + (2^2 + a 2 + /?i) 2 , 



and a new integral G pointed out in [3], which in case (3 = turns into the 
square of the vertical component of the angular momentum. 

Below we exclude the case \a\ = a ■ (3 = 0, for which there exists a cyclic 
coordinate [6], and the problem reduces to the system with two degrees of 
freedom. 

First, we show that without loss of generality one can always take c?_L/3. 
The conditions (2) hold if we replace ei,e 2 ,a,f3 with e*i(6 ) ), e 2 (9), a(8), (3(0), 
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(ei(0),eb(0)) = (?i,? 2 )6, We), (3(9)) = (<2,P)Q, 

, ' cos — sin 
0=1 . ,, „ , = const. 



where 



sin 6 1 cos 

Therefore, e*i(6 l ), e 2 (#) may be taken as the first vectors of a new movable basis. 
At the same time, substitution (4) preserves the rotating moment e*i x a+e 2 x (3 
in Euler equations, and new vectors a(Q),(3(9) satisfy Poisson equations. 

For the general case 

(\a\ 2 - 0\ 2 ) 2 + (a ■ 0) 2 ^ 

take tan2fl = 2(aJ)/(\a\ 2 - \(3\ 2 ) if \a\ ^ and cos 29 = if \a\ = |/3|. 
Then a(9)-L(3(9). Thus, below we consider the natural restrictions (also called 
the geometrical integrals) in the form 

\a\ 2 = a 2 , \(3\ 2 = b 2 , a- (3 = 0. (5) 

The first integral [3] can then be written in a simple way 

G = \{^ 2 a + uf) + - w 3 u; 7 - 6 2 «i - a 2 /3 2 , (6) 

where 

uj a = 2uj\ol\ + 2uj 2 a 2 + uj 3 a 3 , 
up = 2wi/3i + 2u; 2 ^2 + ^3/^3, 

tu 7 = 2u;i(a 2 /3 3 - a 3 /3 2 ) + 2u 2 (a 3 pi - aM + uJ 3 {ai(3 2 - a 2 Pi). 



The conditions (5) make the phase space diffeomorphic to M 6 = R 3 x 50(3). 
In general, the integral manifold Jh,k,g — {H = h,K = k,G = g} C M 6 
consists of 3-dimensional tori bearing quasiperiodic trajectories dense on each 
torus for almost all values of the integral constants. Therefore, a 4-dimensional 
invariant submanifold, on which the induced system has a structure of the 
integrable system with two degrees of freedom, must reside in the set of critical 
points of the global integral mapping J = H x K x G : M 6 — > R 3 . One 
submanifold of this type was found in [2]: 

M 4 = {K = 0} C M 6 . 

The topological structure of the induced system on M 4 was studied in [3]. 

Below we deal with the case [3], which generalizes the so-called marvellous 
motions of the 2nd and 3rd classes of Appelrot [7j. 
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3 New equations for the integral manifolds 



Changing, if necessary, the order of vectors in the movable basis we can assume 
that a > b. Denote 

2 2 , 7,2 2 2 7,2 

p = a + b , r —a — b 
and consider the set N A C M 6 of critical points of the function 

F = (2G — p 2 H) 2 — r 4 K 

belonging to the level F — 0. 

In order to obtain simple formulae and to establish the correspondence with 
[4], introduce new phase variables (i 2 = —1) 



wx — o>i + ioj 2 , w 2 = Wx, w 3 = u 3 , 
Xx — (ax — fa) +i(oi2 + Pi), x 2 = x l} 
yx = (ax + (3 2 ) +i(a 2 - (h), V2 = Vu 
zx = a 3 + i(3 3 , z 2 = zx. 

The equations (1) can be then written as follows 



(7) 



2w[ = —(wxW3 + Zx), 2w' 2 

x[ = —X\W 3 + ZxWx, 

y'l = -yiw 3 + z 2 wx, 

2z[ = xxw 2 - y 2 wx, 
Here stroke stands for d/d(it). 
The conditions (5) take the form 



w 2 w 3 + z 2 , 2w' 3 = y 2 - yx, 



.X o 



x 2 w 3 - z 2 w 2 , 



V 2 = V2W 3 - ZxW 2 , 



2z' 



-x 2 wx + yxw 2 . 



z 2 + xxy 2 = r 2 , z 2 +x 2 y!=r 2 , 



xix 2 + ?/i2/2 + 2^12:2 = 2p 2 



and F = 0, V 6 F = give 
where 



Fx = y/ X l X 2 w 3 



Fx = 0, F 2 = 0, 
1 



y/XxX 2 



(x 2 ZxWx + XxZ 2 W 2 ), 



1 r x 2 



F 2 = -[f(w 2 l +x 1 )~ X ^(w 2 2 + x 2 )}. 
Z Xx x 2 



(9) 



(10) 
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The equations (10) correspond to the system of invariant relations found in 
[I]. This fact, in particular, reveals the topological nature of strictly analytical 
results Moreover, it straightforwardly proves that N 4 is a subset of the 
phase space invariant under the flow of the dynamical system (1). 

It is easy to see that almost everywhere on N 4 the system of equations (10) 
has rank 2, so iV 4 has a natural structure of 4-dimensional manifold except, 
maybe, for a thin subset defined by x\x 2 = 0. Fix the constants h, k, g of the 
integrals (3), (6) and introduce new constants 



(the sign of I is arbitrary). Then from the first integrals, conditions (9) and 
equations (10) we obtain on N 4 



m = — 





x 2 Xi 
Z\W\ z 2 w 2 



w 3 = 1 , 

Xi x 2 



(11) 




) — Ix + Jr 4 m 2 — r 2 m(xi + x 2 ) + x 2 = 0, 





(12) 



The equations (11) of the integral manifold Jh,k,g C iV 4 show that in general 
case for given m, £ this manifold is two-dimensional. 



4 Separation of variables 



We now introduce new variables in (x, z)-plane 



si = 



x 2 + z 2 + r 2 
2x 



s 2 = 



x 2 + z 2 — r 2 



2x 



(13) 



Calculating time derivatives from (8) we obtain 




(14) 
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Let 



*(si, s 2 ) = 4msis 2 - 2£(si + s 2 ) + — (f - 1), 

1 m 

= 4ms 2 - 4£s + — (£ 2 - 1). 
m 



Then, having the obvious identity 

V 2 (s 1 ,s 2 )-<5>(s 1 )<5>(s 2 )=4(s 1 -s 2 ) 2 , 
we find from (9), (11), (12) 



Xl = ~ 2 ( gl -g 2 )2 ^ (si ' S2) + V /$ ( s i) $ ( s 2)]> 

X2 = ~ 2( gl -s 2 ) 2 ^ (Sl ' S2) " V^O^S)]' 
(2 SlS2 - p 2 ) - 2 v /( S 2 -a 2 )( S 2 -6 2 ) 



1/2 = 2 



*(si, s 2 ) - v /$(sQ$(g 2 ) 
( 2siS2 - p2) + 2v /( s 2 - a 2 ) ( S 2 -6 2 ) 



*(si,s 2 ) + J$(si)$(s 2 ) 



^i = ^— (>/«?- « 2 + V^- 62 )' (15) 

si - s 2 v v 

z 2 = -^—^si-ai-^st-V), 
si - s 2 ■ 



= r- 



*(si,s 2 ) - v/$(si)$(s 2 ) 



/$(s 2 ) + v /$(s 1 ) 

w 2 = r- 



*(si,s 2 ) + V$(si)$(s 2 ) 

1 



W3 = ^^hA 3 *-* 2 )^) - v^ 2 -" 2 )^)]- 



Substitution of the latter expressions for Xj,Zj,Wj(j = 1,2) into (14) allows 
to obtain the differential equations for si, s 2 in the real form 



dsi 1 f— ^7~r, — r ds 2 1 



(// 2 y (a 2 - S 2 )$( Sl ), ^ = ^ - *1)*(*2). (16) 

Thus, si, s 2 are easily found as elliptic functions of time. The explicit formulae 
for the phase variables Uj, a.j, (3j(j = 1, 2, 3) immediately follow from (7), (15). 
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5 The types of solutions 



The conditions (9) define the global area for the variables (13): 

|si| > a, \s 2 \ < b. 

So, the domain of oscillations (16) for the fixed values of m, £ is obtained from 
the inequalities 

$(si) < 0, $(s 2 ) > 0. 
Bifurcations of solutions (16) with respect to the parameters m, t take place 
when $(±a) = 0, or $(±6) = 0. It leads to a set of lines in (m, -plane 

t = -2am ± 1, i = 2am ± 1, t = -2bm ± 1, I = 2bm ± 1. 

Analyzing the evolution of roots of the polynomial $(s), we obtain all different 
types of motions. 

The critical motions appear in the cases when one of the variables Si,s 2 re- 
mains constant, coinciding with the double root of polynomial product in the 
right part of the corresponding equation (16). It obviously leads to the motion 
of the body of pendulum type: either u)\ = u 3 = 0, or u 2 = uj 3 = 0. 

An interesting aspect of the case considered is that both movable and immov- 
able hodographs of the angular velocity are explicitly found simultaneously 
without any further integration. Actually, because for almost all constants 
m, £ from the domain of existing of real solutions, hodographs fill some two- 
dimensional surfaces densely, we use the expressions (15) to obtain the para- 
metric equations of those surfaces, in which si, s 2 are independent parameters. 

— * — * 

Expressions for a, (3 (and, therefore, for a x (3) via Si,s 2 give the paramet- 
ric equations for the orientation matrix. Contemporary methods of computer 
graphics provide the possibility to construct a detailed and clear picture of 
motion as rolling without slipping of one surface through the other: at any 
moment these surfaces have a common point with zero absolute velocity and 
one common tangent vector, expressing the fact that absolute and relative 
derivatives of the angular velocity coincide. 
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